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ABSTRACT 

Basement depth estimation of Sokoto Sedimentary Basin was carried out to determine the sedimentary thickness 
in all parts of the basin, The total magnetic field values of the basin, obtained by digitizing the contour maps of 
the basin, which was used to produce the composite map of the area. The Sokoto Basin is situated between 
latitudes 10" 30" N to 14° 00" N and longitude 3 0 30" E to 7° 00" E with an estimated area of 59,570 
km 2 .Spectral depth analysis was carried out on each of the twenty nine aeromagnetic sheets covering the basin, 
to estimate the depths to basement. The results of the spectral studies indicate an increase in sedimentation 
northwards, with several depressions on the basement rock. Two prominent magnetization layers of depths 
varying from 0.04 km to 0.95 km and 0.36 km to 2.88 km were observed. The areas, where higher sedimentary 
thicknesses are observed such as sheets No 31 (TalataMafara), sheet no. 12 (Isa) and sheet no. 13(Gandi) are 
the most probable sites for prospect of hydrocarbon accumulation in the basin. 
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I. INTRODUCTION 

The earth and its contents have long been of concern to mankind. Man has tried to unravel its 
complexity and delve into its origin via various geophysical methods. The subsurface has been of particular 
concern to geoscientists, who seek to investigate it using diverse means, some for the purpose of having 
knowledge, while others do it for exploration of economic resources such as minerals and hydrocarbons. With 
the advances in technology and the need to have a clearer picture of the earth subsurface and its contents, the 
earth scientists have deemed it necessary to utilize the properties associated with earth's interior. Geophysics 
involves the application of physical principles and quantitative physical measurements in order to study the 
earth's interior, its atmosphere and terrestrial space. The analysis of these measurements can reveal howthe 
earth interior varies both vertically and laterally, and the interpretation of which can reveal meaningful 
information on the geological structures beneath ,Dobrin(1976). By working at different scales, geophysical 
methods may be applied to a wide range of investigations from studies of the entire Earth to exploration of a 
localized region of the upper crust for engineering or other purposes (Keareyet al., 2004). A wide range of 
geophysical methods exist for each of which there is an operative physical property to which the method is 
sensitive. The type of physical property to which a method responds clearly determines its range of application. 

Thus, for instance, magnetic method is very suitable for locating buried magnetic ore bodies because 
of their magnetic susceptibility. Similarly, seismic and electrical methods are suitable for locating water table, 
because saturated rock may be distinguished from dry rock by its higher seismic velocity and higher 
electrical conductivity (Keareyet al., 2004). In exploration for subsurface resources, the geophysical methods are 
capable of detecting and delineating local features of potential interest. Geophysical methods for detecting 
discontinuities, faults, joints and other basement structures, include the following: magnetics, seismic, 
resistivity, electrical, potential field, well logging, gravity, radiometric, thermal etc (Corell andGrauch,1985). 
Some geophysical methods such as gamma-ray spectrometry and remote sensing measure surface attributes; 
others, such as thermal and some electrical methods are limited to detecting relatively shallow subsurface 
geological features. Geophysical modeling provides generalized and no-unique solution to questions concerning 
the concerning the geometry of the subsurface geologic structures (Reeves, 2005). Most economic minerals, oil, 
gas, and groundwater lie concealed beneath the earth surface, thus hidden from direct view. 
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The presence and magnitude of these resources can only be ascertained by geophysical investigations 
of the subsurface geologic structures in the area. If the area under investigation has no previous geological 
information and the primary aim of the study is to search for hydrocarbon deposits; the first question that must 
be answered, is whether the sedimentary basin is large enough and thick enough to justify any further 
investigations ( Reynolds, 1990). Other criteria for assessing a basin for possible hydrocarbon accumulation are 
as follows: Source rock: The presence of a source rock with adequate organic materials of various types 
determines whether oil or gas will be sourced in the process of hydrocarbon generation. Reservoir rock: 
Reservoir rock ensures that hydrocarbon generated are accommodated in pore spaces that exists between the 
grains in the subsurface. 

Seal: seals are normally without pore spaces ensures that hydrocarbon generated in the reservoir rocks are kept 
in place and prevented from migrating to the surface and subsequent loss. 

Paleotemperatures: Paleotemperatures or ancient Paleotemperatures are responsible for generation of 
hydrocarbon from organic matter. 

Trap:The trap ensures that oil sealed in place does not migrate and be lost. 

In this study, aeromagnetic data of theSokoto Basin is analyzed by techniques of spectral depth analysis in 
order to estimate the sedimentary thickness in all parts of the basin and suggest areas of likely hydrocarbon 
accumulation prospects.. 

II. MATERIALS AND METHODS 

2.1 Data Acquisition 

The study area (Sokoto Sedimentary Basin) is covered by thirty aeromagnetic maps of total field 
intensity in half -degree sheets. These maps are obtained from the Nigerian Geological Survey Agency 
(NGSA), now renamed the Nigerian Geological Survey Agency (NGSA)).The Agency carried out an airborne 
magnetic survey of substantial part of Nigeria between 1974 and 1980. The magnetic information consists of 
flight lines plotted on a continuous strip chart or tape records. The data were collected at a nominal flight 
altitude of 152.4m along N-S flight lines spaced approximately 2 km apart. The magnetic data collected were 
published in the form Vi degree aeromagnetic maps on a scale of 1:100,000. The magnetic values were plotted at 
lONt (nano Tesla) interval. The maps are numbered, and names of places and coordinates (longitude and 
latitudes) written for easy reference and identification. The actual magnetic values were reduced by 
25,000gamma before plotting the contour map. This implies that the value 25,000 gamma is to be added to the 
co ntour values so as to obtain the actual magnetic field at a given point. A correction based on the International 
Geomagnetic Reference Field, (IGRF,) and epoch date January 1, 1974 was included in all the maps. The visual 
interpolation method, that is the method of digitizing on Grid Layout was used to obtain the data from field 
intensity aeromagnetic maps covering the study area. The data from each digitized map is recorded in a 19 by 19 
coding sheet which contains the longitude, latitude and the name of the town flown and the sheet number. The 
unified composite dataset for the study area was produced after removing the edge effect. Surfer 8 Geosoft ware 
was used to import the dataset. The dataset consists of three columns (longitude,latitude and magnetic values). 
The composite map was produced using Oasis Montaj.. 

2.2 Regional-residual Separation. 

Magnetic data observed in geophysical surveys comprises of the sum of all magnetic fields produced 
by all underground sources. The composite map produced using such data, therefore contains two important 
disturbances, which are different in order of sizes and generally super-imposed. The large features generally 
show up as trends, which continue smoothly over a considerable distance. These trends are known as regional 
trends. Super-imposed on the regional field, but frequently camouflaged by these, is the smaller, local 
disturbances which are secondary in size but primary in importance. These are the residual anomalies. They may 
provide direct evidence of the existence of the reservoir type structures or mineral ore bodies. 

2.3 Production of Regional and Residual Maps 

The residual magnetic field of the study area was produced by subtracting the regional field from the 
total magnetic field using the Polynomial fitting method. The computer program Aerosupermap was used to 
generate the coordinates of the total intensity field data values. This super data file, for all the magnetic values 
was used for production of composite aeromagnetic map of the study area using Oasis Montaj software version 
7.2 A program was used to derive the residual magnetic values by subtracting values of regional field from the 
total magnetic field values to produce the residual magnetic map and the regional map. 
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2.4 Spectral Depth Analysis 

Determination of depths to buried magnetic rocks is among the principal applications of an 
aeromagnetic data. The depths are commonly computed from measurement made on the widths and slopes of an 
individual anomaly of the aeromagnetic profiles. The statistical approach has been found to yield good 
estimates of mean depth to basement underlying a sedimentary basin(Hahnetal., 1976, Udensi,2001). Spector 
(1968) and Spector and Grant (1970) developed a depth determination method which matches two dimensional 
power spectral calculated from gridded total magnetic intensity field data with corresponding spectral obtained 
from a theoretical model. For the purpose of analyzing aeromagnetic data, the ground is assumed to consist of a 
number of independent ensembles of rectangular, vertical sided parallelepiped, and each is ensemble 
characterized by a joint frequency distribution for the depth (h) and length( b) and depth extent ( t).In this work, 
the characteristics of the residual magnetic field are studied using statistical spectral methods. This is done by 
first transforming the data from space to the frequency domain and then analyzing their frequency 
characteristics. In the general case, the radial spectrum may be conveniently approximated by straight line 
segments, the slopes of which relate to depths of the possible layers, (Spector and Grant, 1970, and Hahn 
etal.,1976). The residual total magnetic field intensity values are used to obtain the two dimensional Fourier 
Transform, from which the spectrum is to be extracted from the residual values T (x,y) consisting of M rows and 
N columns in X - Y. The two dimensional Fourier transforms is obtained as given in equation (3.1) above. The 
evaluation is done using an algorithm that is a two dimensional extension of the fast Fourier transform 
(Oppenheim and Schafer, 1975). 

Next, the frequency intervals are subdivided into sub -intervals, which lie within one unit of frequency 
range. The average spectrum of the partial values together constitutes the redial spectrum of the anomalous field 
(Hahn et al., 1976, Kangkolo, 1996 and Udensi, 2001). The logarithm of the energy values versus frequency on a 
linear scale was plotted and the linear segments located; The use of Discrete Fourier Transform introduced the 
problem of aliasing and the truncation effect (or Gibbs phenomenon). Aliasing was reduced by the digitizing 
interval used in the study. Three or two linear segments could be seen from the graphs. The first points on the 
frequency scale was ignored because the low frequency components in the energy spectrum are generated from 
the deepest layers whose locations are most likely in errors, each linear segment groups points due to anomalies 
caused by bodies occurring within a particular depth. If the z is the mean depth of the layer, the depth factor for 
this ensemble of anomalies is exp (-2zk). Thus the logarithmic plot of the radial spectrum would give a straight 
line whose slope is -2z.The mean depth of the burial ensemble is thus given as Where (m) is the slope of the best 
fitting straight line Equation 2.1 can be applied directly if the frequency unit is in radian per kilometer. If 
however, the frequency unit is in circle per kilometer, the corresponding relationship can be expressed as2.2 

III. RESULTS 

3.1 Total Magnetic Intensity Map. 

The total magnetic intensity map (TMI) of the Sokoto sedimentary basin produced from this study 
using Oasis Montaj version 7.2 is as shown in Figure l.The TMI map of Sokoto sedimentary basin can be 
divided into main three sections, though minor depressions exist scattered all over area. The northern part of the 
Sokoto basin is characterized by low magnetic intensity values represented by dark-green-blue color. Whereas 
the southern part is dominated by high magnetic intensity values indicated by red color. The two sections are 
separated by a zone characterized by medium magnetic intensity values area depicted by yellow-orange color. 
These high magnetic intensity values, which dominate the southern part of the Sedimentary basin are caused 
probably by near surface igneous rocks of high values of magnetic susceptibilities. The low amplitudes are most 
likely due to sedimentary rocks and other non-magnetic sources. In general, high magnetic values arise from 
igneous and crystalline basement rocks, whereas low magnetic values are usually from sedimentary rocks or 
altered basement rocks. The sedimentary thickness of theSokoto basin in general, appears to increases from 
south to north. This collaborates well with earlier findings of 2D Seismic surveys conducted by ELF and Mobil 
Companies carried out in 1979. Figure lTotal magnetic intensity map of the study area. 

3.2 Regional Magnetic Intensity Map. 

The regional magnetic intensity map of the study area (Figure 2) produced using Oasis Montaj version 
7.2 is as shown. Figure 2 Regional Map of the Study AreaThe regional magnetic values ranges from 7830 nano 
tesla to 7870 nano tesla and the values decreases from south to north indicating there is a fill of sediments more 
in the northern part of the basin than in the southern part of the study area 
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3.3 Residual Magnetic Intensity Map. 

Figure 3 and Figure4 are the residual magnetic intensity maps of the study area obtained from the total 
magnetic intensity map produced using Surfer 8 and Oasis Montaj version 7.2 respectively. The magnetic 
intensity values ranges from -10 nano tesla to 40 nano tesla. Negative magnetic intensity values are more 
predominant in the northern section of the study area while the southwest has more of positive magnetic 
intensity values. Northeast -Southwest trends are observed in the north central part of the TMI map. 

3.4 Spectral Depth Analysis 

The results for the spectral depths estimates (HI and H2) are as tabulated (Tablel) below. 
The depth estimates obtained are superimposed on the total magnetic intensity map of the 
study area as shown in Figure 5 

TablelTHE DEPTH ESTIMATES 

The contour of basement depth (H2) of the study area is as shown in Figure 6 

The depth contours depicts a number of basement depressions over the entire study area. The depth values 
suggest that the magnetic basement surface beneath the Sokoto sedimentary basin is generally shallow and of 
low relief. Depths to magnetic basement range from 0.36 km in the southern part to 2.85 km at the northwestern 
bordering with Niger republic. A closer observation indicates a general increase in thickness from the southern 
part to northern part, more especially the northwestern part bordering with Niger republic. This general trend in 
sedimentation implies that the Iullemmenden basin is thickest at the centre, which is at Niger republic. 
The surface plot of the second layer (H2) of the study area (Figure 6) shows similar trend with the total 
magnetic intensity map. 

The surface shows the variation in depth of the magnetic source rocks. The peaks shows areas where 
sedimentation are relatively higher. 

Two prominent sedimentary reservoirs are observed in north east and north central part of the map. The result 
agrees well with the results of the upward continuation at depths of 5km and above. 
Figure 6The contour of basement depth (H 
2 

) of the study area. 

DEPTH ESTIMATES OF THE STUDY AREA 

III. CONCLUSIONS 

The conclusion to be drawn from this work, based on the results of spectral depth analysis and 
modeling of the residual magnetic anomalies identified on the Total Magnetic Intensity (TMI) of 
Sokotosedimentary Basin are as follow From this study, results obtained from the spectral depth analysis 
carried out sheet by sheet identified areas of Argungu (2.74km), Isah (2.50km), Mafara(2.67km), Ankah 
(2.50km) Gummi(2.29km),Shanga (2.84km) and Sokoto (2.65km) having re latively higher 
sedimentation. Results of the 2D Modeling carried out on the residual magnetic anomaly map encountered a 
maximum depthof 2.74km in the vicinities of Argungu on profile 2. Several depressions or variations in 
thickness have been observed in some parts of the Sokoto Basin, particularly around Argungu , Isa, Rabah and 
Gandi areas, These deeper sections of the Sokotosedimentary basin identified in this study, might be probable 
potential sites for hydrocarbon deposits and is therefore recommended to be subjected to further investigation. 
Exploration of the Nigerian inland basins is worth given a push. Hydrocarbons if discovered and harnessed 
willincrease the country's reserve and boost productivity. All these will have economic and strategic benefits for 
the country. 
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